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ABSTRACT 

We model gravitational instability in a wide range of isolated disk galaxies, using GADGET, a 
three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter 
halo and a disk of stars and isothermal gas. Absorbing sink particles are used to directly measure the 
mass of gravitationally collapsing gas. Below the density at which they are inserted, the collapsing gas 
is fully resolved. We make the assumption that stars and molecular gas form within the sink particle 
once it is created, and that the star formation rate is the gravitational collapse rate times a constant 
efficiency factor. In our models, the derived star formation rate declines exponentially with time, 
and radial profiles of atomic and molecular gas and star formation rate reproduce observed behavior. 
We derive from our models and discuss both the global and local Schmidt laws for star formation: 
power-law relations between surface densities of gas and star formation rate. The global Schmidt law 
observed in disk galaxies is quantitatively reproduced by our models. We find that the surface density 
of star formation rate directly correlates with the strength of local gravitational instability. The local 
Schmidt laws of individual galaxies in our models show clear evidence of star formation thresholds. 
The variations in both the slope and the normalization of the local Schmidt laws cover the observed 
range. The averaged values agree well with the observed average, and with the global law. Our 
results suggest that the non-linear development of gravitational instability determines the local and 
global Schmidt laws, and the star formation thresholds. We derive from our models the quantitative 
dependence of the global star formation efficiency on the initial gravitational instability of galaxies. 
The more unstable a galaxy is, the quicker and more efficiently its gas collapses gravitationally and 
forms stars. 

Subject headings: galaxy: evolution — galaxy: spiral — galaxy: kinematics and dynamics — galaxy: 
ISM — galaxy: star clusters — stars: formation 



1. INTRODUCTION 

Stars form at widely va rying rates in different disk 
galaxies l)KennicuttJ fl99 8a'l . However, they appear to 
follow two simple empirical laws. The first is the correla- 
tion between the star formation rate (SFR) density and 
the gas density , the "Schmidt law" as first introduced by 
IScnmidtl ll795fr: 



J SFR 



= A E£ 



(l) 



where Esfr and E gas are the surface densities of SFR 
and gas, respectively. 

When E gas and Esfr are averaged over the entire star 
forming regio n of a galaxy, th ey give rise to a global 
Schmidt law. iKennicutH (|199 8bl found a universal global 
star formation law in a large sample that includes 61 
normal spiral galaxies that have Hq,, Hi and CO mea- 
surements and 36 infrared-selected star burst galaxies. 
The observations show that both the slope N ~ 1.3- 
1.5 and the normalization A appear to be remarkably 
consistent from galaxy to gala xy. There are s ome v ari- 
ations, though. For example, iWong fc Blitzl (|2002|) re- 
ported N ~ 1.1-1.7 for a sample of seven molecule-rich 
spiral galaxies, depending on the correction of the ob- 
served H Q emi ssion for extin c tion i n deriving the star 
formation rate. iBoissier et al.l 1)2003?) examined 16 spiral 
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galaxies with published abundanc e gradients and found 
N ~ 2.0. IGao fc Solomonl l)2004aTl surveyed HCN lumi- 
nosity, a tracer of dense molecular gas, from 65 infrared 
or CO-bright galaxies including nearby normal spiral 
galaxies, luminous infrared galaxies, a nd ultraluminous 
infrared galaxies. Based on this survey. IGao fc Solornonl 
l)2004b|) suggested a shallower star formation law with a 
power-law index of 1.0 in terms of dense molecular gas 
content. 

When E gas and Esfr are measured radially within 
a galaxy, a local Schmidt law can be measured. 
Wong & Blitz ( 2002) investigated the local Schmidt laws 
of individual galaxies in their sample. They found simi- 
lar correlations in these galaxies but the normalizations 
and slopes vary from galaxy to galaxy, with N ~ 1.2-2.1 
for total gas, assuming that extinction depends on gas 
column de nsity (or, N ~ 0.8-1 .4 if extinction is assumed 
constant). iHever et al.l l)2004)) reported that M33 has a 
much deeper slope, N ~ 3.3. 

The second empirical law is the star formation thresh- 
old. Stars are observed to f orm efficiently only above a 
critical gas surface density iMartin fc Kennicuttl 1)200 If ) 
studied a sample of 32 nearby spiral galaxies with well- 
measured H a and H2 profiles, and demonstrated clear 
surface-density thresholds in the star formation laws in 
these galaxies. They found that the threshold gas den- 
<sM$ p leasured at the outer threshold radius where SFR 
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drops sharply) ranges from 0.7 to 40 M Q pc~ 2 among 
spiral galaxies, and the threshold density for molecu- 
lar gas is ~5-10 M Q pc -2 . However, they found that 
the ratio of gas sur face den si ty at the threshold to the 
critical density for iToomrel l|1964ft gravitational insta- 
bility, ctQ = Egas/I] C iit, is remarkably uniform with 
oiq = 0.69 ± 0.2. They assumed a constant veloc- 
ity dispersion of the gas, the effective sound speed, of 
c s = 6 km s _1 . Such a density threshold app lies to 
normal disk galaxie s fe.g.. | B oissier et all 12003ft . ellip- 
tical galaxies (e.g., iVader fe Vigroux I |1991l) 7low sur- 
face brightness gal axies llvan der Hulst et al. Ill993() , and 
starburst galaxies l|Elmegreenlll994ft . However, there are 
a few exceptions, such as dwarf and ir regular galaxies 
(e.g., iHunter. Elmegreen fe Baker] IT9 98). Furthermore, 
inefficient star fo rmation can be found well outside the 
threshold radius (Ferguson ei~aIIIT99l . 

What is the origin of the Schmidt laws and the star 
formation thresholds? The mechanisms that control star 
formation in galaxies, such as gravitational instability, 
supersonic turbulence, ma gnetic fields, and rotational 
shear are widely debated IShu et alJ fl 987t lElmegreenl 
120021 lLarsonll2003i iMac Low fc Klessenll2004ft . At least 
four types of models are currently discussed. The first 
type empha sizes sel f-gravi t y of the g alactic disk (e.g., 
i Quirkl 119721: iLarsonl 119 881 iKennicutt) fl98Ht lElmegreenl 
119941 IKennicutt) 11998b 1 ). In these models, the Schmidt 
laws do not depend on the local star formation process, 
but are simply the results of global gravitational collapse 
on a free-fall time. In the second type, the global star 
formation rate scales with either the loc al dynamical 
time, invoking clo u d-cloud co l lisions (e.g . iWvsd 119861 
IWvsefc Silklfl987t HM It^Eopffi. or the local 
orbital time of the galactic disk (e.g. lEhnegreeri 119971 
IHunter et alJ ll9981. A third type, which invokes hi- 
erarchical star for mation triggered b y turbulence, has 
been proposed by Fl megreenl l)2002ft . In this model, 
the Schmidt law is scale-free, and the star formation 
rate depends on the probability distribution function 
(PDF) of the gas density produced by galactic turbu- 
lence, which appears to be log-normal in simulations 
of tur bulent molecular cl o uds and interstellar medium 
(e.g., iScalo et al. I 1199 81 iPasso t fc Vaz auez-Sem acjeni) 
iaal_ [Ostriker. Oammie. fc Stond 
Klessenl 120001 iWada fe Norman I 200 



Ballestcros-Parcdcs & Mac Low J [2002; 

Padoan fc Nordlundl 120021 ILL Klessen fc Mac Low! 
20031 lKravtsovll2003t [Mac Low et al. 112005ft . Recently. 
Krumholz & McKec (2005) extended this analysis with 



additional assumptions such as the virialization of the 
molecular clouds and star formation efficiency to derive 
the star formation rate from the gas density PDFs. 
They successfully fitted the global Schmidt law, but 
their theory still contains several free or poorly con- 
strained parameters, and does not address the observed 
variation in local Schmidt laws among galaxies. A fourth 
type appeals to the gas dynamics and thermal state 
of the gas to determ ine the star formation behavior. 
iStruck-Marcelll l)1991ft and citetstruck99, for example, 
suggest that galactic disks are in thermohydrodynamic 
equilibrium maintained by feedback from star formation 
and countercirculating radial gas flows of warm and cold 
gas. 

There is considerable debate on the star formation 



threshold as well. iMartin fc Kennicuttl l)2001ft suggest 
that the t hreshold densit y is determined by the Toomre 
criterion (IToomrel 11964^) for gravitational instability. 
IHunter et al.1 (|199q) argued that the critical density for 
star formation in dwarf galaxies depen d s on the rota- 
tional shear of the disk. I Wong fc Blitz) l)2002ft claimed 
no clear evidence for a link between a.Q and star for- 
mation. Instead, they suggested tha t ap is a measure- 
ment of gas fraction. iBoissier et alJ l)2003ft found that 
the gravitational instability criterion has limited appli- 
cation to their sample. Note all the models above are 
base d on an assum ption of constant sound speed for the 
gas. iSchavel Ip0 041 proposed a thermal instability model 
for the threshold, in which the velocity dispersion or ef- 
fective temperature of the gas is not constant, but drops 
from a warm (i.e. 10 4 K) to a cold phase (below 10 3 K) 
at the threshold. He suggested that such a transition is 
able to reproduce the observed threshold density. 

While each of these models has more or less succeeded 
in explaining the Schmidt laws or the star formation 
threshold, a more complete picture of star formation on a 
galactic scale remains needed. Meanwhile, observations 
of other properties related to star formation in galaxies 
have provided more clues to the dominant mechanism 
that controls global star formation. 

An analysis of the distribution of dust in a 
sample of 89 edge-on, bulgeless d isk galaxies by 
iDalcanton. Yoachim fc Bernstein! l)2004ft shows that dust 
lanes are a generic feature of massive disks with V ro t > 
120 km s -1 , but are absent in more slowly rotating 
galaxies with lower mass. These authors identify the 
V I0 t — 120 km s _1 transition with the onset of gravi- 
tational instability in these galaxies, and suggest a link 
between the disk instability and the formation of the dust 
lanes which trace star formation. 

Color gradients in galaxies help trace their star forma- 
tion history by revealing the distribution of the ir stellar 
populations ijSearle. Sargent fc Bagnuolo Il973j) . A com- 
prehensive study of color g r adien ts in 121 nearby disk 
galaxies by iBell fc de .Tonal l)2000ft shows that the star 
formation history of a galaxy is strongly correlated with 
th e surface mass dens i ty. Si milar conclusions were drawn 
by Kauffmann et al. (2003) from a sample of over 10 5 
galaxies from the Sloan Digital Sky Survey. Recently, 
iMacArthur et all l)2004ft carried out' a survey of 172 low- 
inclination galaxies spanning Hubble types SO-Irr to in- 
vestigate optical and near-IR color gradients. These au- 
thors find strong correlations in age and metallicity with 
Hubble type, rotational velocity, total magnitude, and 
central surface brightness. Their results show that early 
type, fast rotating, luminous, or high surface brightness 
galaxies appear to be older and more metal-rich than 
their late type, slow rotating, or low surface brightness 
counterparts, suggesting an early and more rapid star 
formation history for the early type galaxies. 

These observations show that star formation in disks 
correlates well with the properties of the galaxies such as 
rotational velocity, velocity dispersion, and gas mass, all 
of which directly determine the gravitational instability 
of the galactic disk. This suggests that, on a galactic 
scale, gravitational instability controls star formation. 

The nonlinear development of gravitational instability 
and its effect on star formation on a galactic scale can be 
better understood through numerical modeling. There 
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have been many simulations of disk galaxies, including 
isolated galaxies with various a ssumptions of the gas 
physi c s and feedback effects (e. g ..iThacker fc CouchmarJ 
200(1 IWada fc Normanl l200lt INoguchil 120011: IBarnesj 
IRobertson et "aTT"l2004t ILi. Mac Low fc Klessenl 
lOkamo to et al l 120051). galaxy mergers (e.g., 
Mihos ~ Hernauistl 119941 iBarnes fc Hernauistl Il996t 



Li. Mac Low fe Klessenl |2004[) . and galaxies in a cos- 



mological context, with different assumptions about 
the nature and distribution of 




dark 
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matter 
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e.g., 
KatJ 1992 



Katz fc Gu nn 19 911 INavarro fc Bend 119911 

inmetz fc Mueller! [19941 INavarro. Frenk fc Whit 
ISommer-Larsen. Gelato fc Veda 
iSteinmetz fc Navarro! 11 999: iSpr inge 
iSom rner^Larsen fc Dolgovl [200 

Sommer-Larsen.^5t^^^PorthTaril E 

Springel fc Hernauistl 120031 iGovernato et al.l [200 
However, in most of these simulations gravitational 
collapse and star formation are either not numerically 
resolved, or are followed with empirical recipes tuned to 
reproduce the observations a priori. There are only a 
handful of numerical studies that focus on the star for- 
mation laws. Early three-dimensional smoothed particle 
hydrodynamics (SPH) simulati ons of isolated barre d 
galaxies were carried ou t bv IFriedli fc Benzl f l993~): 
IFriedli. Benz fc Kenniaittl ill 9941) . and iFriedli fc Benz I 
(1995). In IFriedli fc Bena l[1993fl . the secular evolution 
of the isolated galaxies was followed by modeling of 
a two-component (gas and stars) fluid, restricting the 
interaction between the two to purely gravitational 
coupling. Th e simu lations were improved later in 
IFriedli fc Benz I (^995) by including star formation 
and radiative cooling. These authors found that their 
method to simulate star formation, based on Toomre's 
criterion, naturally reproduces both the density thresh- 
old of 7 M Q pc~ 2 for star formation, and the global 
Schmidt law in disk galaxies. They also found that the 
nuclear starburst is associated with bar for mation in 
the galactic center. iGerritsen fc Ickd lll997j) included 
stellar feedback in similar two-component (gas and 
stars) simulations that yielded a Schmidt law with 
power-law index of ~ 1.3. However, these simulations 
included only star s and gas, and no dark matter. More 
recently. iKravtsovl <[2003l) reproduced the global Schmidt 
law using self-consistent cosmological simulations of 
high-redshift galaxy formation. He argued that the 
global Schmidt law is a manifestation of the overall 
density distribution of the interstellar medium, and 
that the global star formation rate is determined by 
the supersonic turbulence driven by gravitational in- 
stabilities on large scales, with little contribution from 
stellar feedback. However, the strength of gravitational 
instability was not directly measured in this important 
work, so a direct connection could not be made between 
instability and the Schmidt laws. 

In order to investigate gravitational instability in disk 
galaxies and consequent star formation, we model iso- 
lated ga laxies with a wid e range of masses and gas frac- 
tion. In lLi et aD l[2005al hereafter Paper I) we have de- 
scribed the galaxy models and computational methods, 
and discussed the star formation morphology associated 
with gravitational instability. In that paper it was shown 
that the nonlinear development of gravitational insta- 
bility determines where and when star formation takes 



place, and that the star formation timescale tsf depends 
exponentially on the initial Toomre instability parame- 
ter for the combination of colliso nless stars and collisional 
gas in the disk Q sg derived bv iRafikovl ((2001). Galaxies 
with high initial mass or gas fraction have small Q sg and 
are more unstable, forming stars quickly, while stable 
galaxies with Q sg > 1 maintain quiescent star formation 
over a long time. 

Paper I emphasized that to form a stable disk and 
derive the correct SFR from a numerical model, the 
gravitational collapse of the gas must be fully resolved 
(Bate fc Burkerd 119971: iTruelove et~aTlll997[) up to the 
density where gravitationally collapsing gas decouples 
from the flow. If this is done, stable disks with SFRs 
comparable to observed values can be derived from mod- 
els using an isothermal equation of state. With insuffi- 
cient resolution, however, the disk tends to collapse to the 
center producing much higher SFRs, as f ound by some 
previous work (e.g. IRobertson et alJl2004lh 

We analyze the relation between the SFR and the gas 
density, both globally and locally. In §0wc briefly review 
our computational method, galaxy models and parame- 
ters. In §13 we present the evolution of the star formation 
rate and radial distributions of both gas and star forma- 
tion. We derive the global Schmidt law in § 0[ followed 
by a parameter study and an exploration of alternative 
forms of the star formation law. Local Schmidt laws are 
presented in §[SJ In §[S]we investigate the star formation 
efficiency. The assumptions and limitations of the mod- 
els are discussed in § Finally, we summarize our work 
in § [S| Preliminary results on the global Schmidt law 
and star formation thresholds w ere already presented by 
ILi. Mac Low fc Klessenl J2005rl. 

2. COMPUTATIONAL METHOD 

We here summarize the algorithms, galaxy mod- 
els, and numerical parameters described in detail 
in Paper I. We use the SPH code GADGET, vl.l 
([Springel. Yoshida fc White) 120011). modified to inclu de 
absorbing sink particles l[Bate. Bonnell fc Pricelll995[) to 
directly measure th e mass of grav i tation ally collapsing 
gas. Paper I and Uappsen et a l. (2005) give detailed 
descriptions of sink particle implementation and inter- 
pretation. In short, a sink particle is created from the 
gravitationally bound region at the stagnation point of a 
converging flow where number density exceeds values of 
n = 10 cm~ 3 . It interacts gravitationally and inherits 
the mass, and linear and angular momentum of the gas. 
It accretes surrounding gas particles that pass within its 
accretion radius and are gravitationally bound. 

Regions where sink particles form have pressures 
P/k ~ 10 7 K cm" 3 typical of massive star-forming re- 
gions. We interpret the formation of sink particles as 
representing the formation of molecular gas and stellar 
clusters. Note that the only regions that reach these high 
pressures in our simulations are dynamically collapsing. 
The measured mass of the collapsing gas is insensitive to 
the value of the cutoff-density This is not an important 
free parameter, unlike in the mo dels of lElmegreenl ( 2002) 
and lKrumholz fc McKeel l[2005[) . 

Our galaxy model consists of a dark matter halo, 
and a disk of stars and isothermal gas. The ini- 
tial __gahixv_stnic_ture isbased on the analytical work 
bv IMo. Mao fc White] (1998J), as implemented numeri- 
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callv bv ISpringel fc White! Jl99<l: iSpringell l|2000(l and 
ISpringel. Di Matteo fe H crnauist (2005). We character- 
ize our models by the rotational velocity V200 a t the virial 
radius -R200 where the density reaches 200 times the cos- 
mic average. We have run models of galaxies with rota- 
tional velocity V200 = 50-220 km s _1 , with gas fractions 
of 20-90% of the disk mass for each velocity. 

Observations of Hi in many spiral galaxies suggest 
that the gas velocity dispersion has a range of ~4- 
15 k m s _1 (e.g., see review in I Dib. Bell fc Burkerd 
2005). The dispersion varies radially from ~ 12- 
15 km s _1 in their cen tral regions to ~4-6 km s -1 in 
the outer parts (e.g.. Ivan der Kruit fc Shostakl IT982t 
Dickev. Hanson, fc Heloul I199C iKamphuis fc S ancisi 
19931 iRownd. Dickev. fc Heloul Il994t iMeurer et al. 
1996). We therefore choose two sets of effective sound 
speeds for the gas, c s = 6 km s -1 (low temperature 
models) as suggested bv lKennicutd l)1998b[K and c s = 15 
km s" 1 (high temperature models). Table 1 lists 
the most important model parameters. The Toomre 
criterion for gravitational instabilit y that couples stars 
and gas, Q sg is calculated following Rafikov (2001), and 
the minimum value is derived using the wavenumber k 
of greatest instability and lowest Q sg at each radius. 

Models of gravitational collapse must satisfy three 
numerical criter i a: t he Jeans resolutio n criterion 
(IBate fc BurkerU 119971 hereafter BB97; IWhitwortbl 
1998), the gravity- hydro balance criterion for gravita- 
tional softening (BB97) , and the equipartition criterion 
for particle masses ijSteinmetz fc White! 11997^ . We set 
up our simulations to satisfy the above three numerical 
criteria, with the computational parameters listed in Ta- 
ble 1. We choose the particle number for each model such 
that they not only satisfy the criteria, but also so that all 
runs have at least 10 6 total particles. The gas, halo and 
stellar disk particles are distributed with number ratio 
N g : A/h : -ZVd = 5 : 3 : 2. The gravitational softening 
lengths of the halo = 0.4 kpc and disk — 0.1 kpc, 
while that of the gas h g is given in Table 1 for each model. 
The minimum spatial and mass resolutions in the gas are 
given by h g and twice the kernel mass (~ 80m 9 ). (Note 
that we use h g here to denote the gravitational soften- 
ing length instead of e as used in previous papers, to 
distinguish it from the star formation efficiency used in 
later sections.) We adopt typical values for the halo con- 
centration parameter c = 5, spin parameter A = 0.05, 
and H ubble constant Ho — 70 km s" 1 Mpc" 1 ijSpringell 
2000). 

Resolution of the Jeans lengt h is vital for simu- 
lations of gravitation al collapse l|Truelove et alJ Il997t 
IBate fc Burkerti 11997). Exactly how well the Jeans 
length mu st be resolved r e mains a point of controversy, 
however. iTruelove et al.l l)1997|) suggest that a Jeans 
mass must be resolved with far more than the Nk = 2 
smoothing kernels proposed by BB97. In this work, we 
carried out a resolution study on low-T models G100-1 
and G220-1, with three resolution levels having total par- 
ticle numbers N tot = 10 5 (Rl), 8 x 10 5 (R8) and 6.4 x 10 6 
(R64), reaching Nk « 24. The particle numbers are cho- 
sen such that the maximum spatial resolution increases 
by a factor of two between each pair of runs. Paper I 
finds convergence to within 10% of the global amount of 
mass accreted by sink particles between the two highest 
resolutions, suggesting that the BB97 criterion is suffi- 



cient for the problem of global collapse in galactic disks. 
In this paper we also refer to the resolution study as 
applicable. 

3. STAR FORMATION AND GAS DISTRIBUTION 

Molecular clouds and stars form together in galaxies. 
Molecular hydrogen forms on dust gr ains in a time of 
l|Hollenbach. Werner, fc Salpeterlll97l|) : 

*/ * (10 9 yr)/n (2) 

where n is the number density of the gas. The absence 
of > 10 Myr old stars in star-forming regions in the 
solar neighborhood suggests that molecular cloud 
complexes must coalesce rapidly and form stars quickl y 
( Ball ester os-Pare des. Hart mann. fc Vazau ez-Semadenil 
ll999HHartmannll2000|) . lHartmann et alJ l)2001|) further 
suggested that the conditions needed for molecular 
gas formation from atomic flows are similar to the 
conditions needed for gravitational instability. Star 
formation can therefore take place within a free-fall 
timescale once molecular clouds a re produced. Using 
a one-dimensional chemical model, iBergin et al.l (2004) 
showed rapid formation of molecular gas in 12-20 Myr in 
sho ck-compressed regions. Th ese results are confirmed 
by lOlover fc Mac Low! l)2005l) using three-dimensional 
magnetohydrodynamics simulations with chemistry of 
supersonic turbulence. They show that most of the 
atomic gas turns into H2 in just a few megayears once 
the average gas density rises above ~ 100 cm" 3 . This 
is because the gas passes through turbulent density 
fluctuations of higher density where H2 can form quickly. 
By the time gas reaches the densities of n = 1000 cm" 3 
where we replace it with sink particles, the molecular 
hydrogen formation timescale tf < 1 Myr (eq. |2J. 
Motivated by these results, we identify the high-density 
regions formed by gravitational instability as giant 
molecular cloud complexes and replace these regions by 
accreting sink particles. 

We assume that a fraction of the molecu lar gas turns 
into s tars quickly. CO observat ions by lYoung et alJ 
(1996) and lRownd fc Yound l)1999|) suggest that the local 
star formation efficiency (SFE) in molecular clouds re- 
mains roughly constant. To quantify the SFR, we there- 
fore assume that individual sink particles form stars at a 
fi xed local efficiency q. 

IKennicutd ll1998hfl found a global SFE of e„ = 30% 
for starburst galaxies, which IWong fc Blitl l)2002j) found 
to be dominated by molecular gas. We take this to be 
a measure of the local SFE ei in individual molecular 
clouds, since most gas in these galaxies has already be- 
come molecular. In our simulations, sink particles rep- 
resent high pressure (P/k « 10 7 K cm" 3 ), massive star 
formation regions in galaxies, such as 30 Doradus i n the 
Large Magellanic Cloud (e.g.. iWalborn et al]l!999|) . We 
therefore adopt a fixed local SFE of q = 30% to convert 
the mass of sink particles to stars, while making the sim- 
ple approximation that the remaining 70% of the sink 
particle mass remains in molecular form. This approach 
will be discussed in more detail in § El 

3.1. Evolution of Star Formation Rate 

FigureQ]shows the time evolution of the SFRs of differ- 
ent models. The SFR is calculated as SFR = dM*/dt ~ 
q A M s i n k/At, where M* and M s i„k are the masses of 
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Fig. 1. — (a) Time evolution of the SFRs in selected high-T 
models as given in the legend, (b) Resolution study of low-T models 
G100-1 (blue) and G220-1 (orange) with resolutions of Rl (dotted 
lines), R8 (dashed lines) and R64 (solid colored lines), where the 
resolution levels are in units of 10 total particles. The standard 
R10 models (black) are shown for comparison. 
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Fig. 2. — Correlation between the maximum of the SFR and the 
minimum radial values of the Toomre instability parameters (a) 
for stars and gas together Q S g jm i n and (b) for gas alone Q g m j n . 
The solid lines are the least absolute deviation fits log SFR max = 
K + AQ min . 



the stars and sink particles, respectively. We choose 
€i = 30% to be the local SFE within sink particles, and 
the time interval At — 50 Myr. Figure^, shows the SFR 
curves of selected high T models, while Figure^) shows a 
resolution study of SFR evolution. We find convergence 
to within 10% of the SFR between the two highest reso- 
lutions over periods of more than 2 Gyr, suggesting that 
this result converges well under the BB97 criterion. 

The star formation rates in Figure Q decline over time. 
Many of the models have SFR > 10M Q yr _1 , cor- 



responding to starburst galaxies. Some small or gas 
poor models such as G100-1, on the other hand, have 
SFR < IMq yr _1 . They maintain slow but steady star 
formation over a long time and may represent quiescent 
normal galaxies. The maximum SFR appears to depend 
quantitatively on the initial instability of the disk as mea- 
sured either by Q S g,min, the minimum Toomre parame- 
ter for the combination of stars and gas in the disk, or 
by the value for the gas only Q g>m i n . Figure |21 shows 
both correlations: logSFR max ~ 3.32 — 3.13Q S9 m i n and 
log SFR max ~ 2.76 - 1.99Q 9 , min . 



() 



The SFRs in most models shown in Figure ^ appear 
to decline exponentially. From Paper I, the accumulated 
mass of the stars formed in each galaxy can be fitted with 
an exponential function: 



M* 



M [1 - exp (-*/tsf)] 



(3) 



where 



Mo = 0.96 q {1 - 2.9 exp [-1.7/Q sa , roin ]} , (4) 

r SF = (34 ± 7 Myr) exp [Q sg , min /0.24] , (5) 

M init is the initial total gas mass, and tsf is the star 
formation timescale. The star formation rate can then 
be rewritten in the following form: 

' (6) 



dM* 

SFR = — - tx exp -t/rsF 

at t sf 



A similar exponential form is also repor ted by 
IMacArthur et al.1 II2004D. as first su ggested by ILarsonl 
(I1974D andlTinslev fc Larson I l(l97^ . 

Sandage] l)1986f) studied the star formation rate of dif- 
ferent types of galaxies in the Local Group and proposed 
an alternative form fo r the star formation hist ory as ex- 
plicitly formulated bv IMacArthur et alJ l|2004|) : 

t 



SFRcx 



exp i 



-t 2 /rl F ) 



(7) 



'SF 



Figure shows an example of model G220-1 (low- 
T) fitted with these two formula. Both formulae appear 
quite similar at intermediate times. The Sandage model 
captures the initial rise in star formation better, but the 
exponential form follows the late time behavior of our 
models more closely. As we only include stellar feedback 
implicitly by maintaining constant gas sound speed, we 
must be somewhat cautious about our interpretation of 
the late time results. In order to compare the fits, we 
define a parameter for relative goodness of the fit 

x 2 = E([^-y/]/^) 2 (8) 

where y s is the SFR from the simulation, y m is the maxi- 
mum of SFR, and yj is the model function from equation 
(JBJ or (J2J. Note that since we do not take into account 
the uncertainty of each point, the absolute value of \ 2 
has no meaning. We only compare the relative \ 2 val- 
ues in Figure Both formulae fit equally well to many 
models, especially to those with high gas fractions that 
form a lot of stars early on. But for some models such 
as G100-1 (low-T) and G220-1 (high-T), the exponential 
function seems to fit noticeably better. Therefore we use 
the exponential form in the rest of the paper. 

This analysis implies that the star formation history 
depends quantitatively on the initial gravitational in- 
stability of a galaxy after its formation or any major 
perturbation. An unstable galaxy forms stars rapidly 
in an early time, so its stellar populations will appear 
older than those in a more stable galaxy. More mas- 
sive galaxies are less stable than small galaxies with the 
same gas fraction. The different star formation histo- 
ries in such galaxies may account for the downsizing ef- 
fect that star formation first occurs in big galaxies at 
high redshift. while modern starburs t gala xies are smaller 
llCowie et al.lll996t IPoggianti et alJl2004t iFerreras et all 
I2004D . and thus more stable. 
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Fig. 3. — (a) Example fit of t he SFR evolutio n curve with the 
exponential form (eq. |SJ and the Sandage ( 1986;) form (eq. for 
model G220-1 (low-T) as an example, (b) The relative goodness 
of fit y 2 (eq. [S) for all models for exponential (red) and ISandage I 
11986) (black) forms. 



3.2. Radial Distribution of Gas and Star Formation 

Figure 01 shows the radial distribution of different gas 
components and the SFR of selected models. The gas 
distribution and SFR are calculated at the star formation 
timescale tsf derived from the fits given in Paper I. We 
assume that 70% of the gravitationally collapsed, high- 
density gas (as identified by sink particles) is in molecular 
form. Similarly, we identify unaccreted gas as being in 
atomic form. The total amount of gas is the sum of both 
components. 

Figure 0] shows that the simulated disks have gas dis- 
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Fig. 4. — Radial profiles of atomic (SPH particles; green dimond), molecular (70% of sink particle mass; blue square), and total (red 
dot) gas surface density, as well as SFR surface density (black dot). R d is the radial disk scale length as given in Table 1. Note the solid 
lines are used to connect the symbols. 



tributions that are mostly atomic in the outer disk but 
dominated by molecular gas in t he ce ntral region. Obser- 
vations bv lWong k Blitil|200^ and lHever et~al] pool 
show that the density difference between the atomic and 
molecular components in the central region depends on 
the size and gas fraction of the galaxy. For example, Sh 2 
in the center of NGC432 1 is almost two orde rs of mag- 
nitude higher than E Hl l|Wong fc Blitz! 12002) , while in 
M33, the difference is only about one order of magnitude 
ijHever et al.ll200 4'). A similar relation between the frac- 
tion of molecular gas and the gravitational instability of 
the galaxy is seen in our simulations. In our most un- 
stable galaxies such as G220-4, the central molecular gas 
surface density exceeds the atomic gas surface density by 
more than two orders of magnitude, while in a more sta- 
ble model like G220-1 (high T), the profiles of S H2 and 
Shi are close to each other within one disk scale length 
i?d- 

We also find a linear correlation between the molecular 



gas surface density and the SFR surface density, as can 
be seen by their parallel radial profiles in Figure (In 
operational terms, we find a correlation between the sur- 
face density in s ink particles and the rat e at which they 
accrete mass. ) IGao fc Solomon! (|2004bD found a tight 
linear correlation between the far infrared luminosity, a 
tracer of the star formation rate, and HCN luminosity, in 
agreement with our result that star formation rate and 
molecular gas have similar surface density profiles. 

The agreement between the simulations and observa- 
tions supports our assumption that both molecular gas 
and stars form by the gravitational collapse of high den- 
sity gas. Note, however, that we neglect recycling of gas 
from molecular clouds back into the warm atomic and 
dissociated or ionized medium represented by SPH par- 
ticles in our simulation. Although it is possible that even 
that reionized gas may still quickly collapse again if the 
entire region is gravitationally unstable, this still con- 
stitutes an important limitation of our models that will 
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have to be addressed in future work. 

4. GLOBAL SCHMIDT LAW 

To derive the global Schmidt law from our models, we 
average Ssfr and S gas over the entire star forming re- 
gion of e ach galaxy. We d efine the star forming region 
following IKennicuttl <|1989j) . using a radius chosen to en- 
circle 80% of the mass accumulated in sink particles (de- 
noted i?80 hereafter). The SFR is taken from the SFR 
evolution curves at some chosen time. As mentioned in 
§ 01 30% of the mass of sink particles is assumed to be 
stars, while the remaining 70% of the sink particle mass 
remains in molecular form. The atomic gas component 
is computed from the SPH particles not participating 
in localized gravitational collapse, that is, gas particles 
not accreted onto sink particles. The total is the sum of 
atomic and molecular gas. 

Figure [S] shows the global Schmidt laws derived from 
our simulations at the star formation time t = tsf as 
listed in Table 1. Note that for a few models this is just 
the maximum simulated time, as indicated in Table 1. 
(Results from different times and star formation regions 
are shown in the next section.) We fit the data to the 
total gas surface density of the models listed in Table 1 
(both low T and high T) . A least-square fit to the models 
we have run gives a simulated global Schmidt law 



-'SFR 



(1.1 ±0.4 x 10" 



M yr 1 kpc 

t 1.56±0.09 

Jgas 



1 M pc" 



(9) 



For comparison , the best fit to the observations by 
IKennicuttl l|1998bf) gives a global Schmidt law for the 
total gas surface density in a sample that includes both 
the normal and starburst galaxies of 



S SF r = (2.5 ± 0.7 x 10~ 4 Mq y^ 1 kpc" 2 ) 

1.4±0.15 



-gas 



1 Mq pC" 



(10) 



The global Schmidt law derived from the simulations 
agrees with the observed slope within the observational 
errors, but has a normalization a bit lower than the ob- 
served range. There are three potential explanations for 
this discrepancy. First, we have not weighted the fit by 
the actual distribution of galaxies in mass and gas frac- 
tion. Second, as we discuss in the next subsection, we 
have not used models at different times in their lives 
weighted by the distribution of lifetimes currently ob- 
served. Third, we have only simulated isolated, normal 
galaxies. Our models there fore do not populat e the high- 
est Ssfr values observed in lKennicutti l|1998b(l . which are 
all starbursts occurring in interacting galaxies. These 
produce highly unstable disks that undergo vigorous 
starbursts with high SFR (e.g.. lLi. Mac Low k, Kless"enl 
2004). Our result is supported bv lBoissier et all {20031. 
who found a much deeper slope, N ~ 2 in a sample of 
normal galaxies comparable to our more stable models. 
In the models, the local SFE q is fixed at 30%, inde- 
pendent of the galaxy model. A change of the assumed 
value of eg changes the normalization but not the slope 
of our relation. For example, an extremely high value of 
ee ~ 90% increases A to ~ 3.15 x 10~ 4 M yr" 1 kpc -2 , 
which is just within the la upper limit of the observation 



bv IKennicuttl l|1998rj|) . If we decrease ee to 10%, then A 
decreases to ~ 0.7 x 10~ 4 M yr -1 kpc -2 . These fairly 
extreme assumptions still produce results lyi ng within 
the observed ranges re.g,. IWong k Blitdl2002V suggest- 
ing that our overall results are insensitive to the exact 
value of the local SFE that we assume. 

The SFR surface densities Ssfr change dramatically 
with the gas fraction in the disk. The most gas-rich mod- 
els (M-4, circles) have the highest Ssfr, while the mod- 
els poorer in gas (M-l, squares) have Ssfr two orders of 
magnitudes lower than their gas-rich counterparts. Note 
that models with lower Ssfr tend to have slightly higher 
scatter, because in these models fewer sink particles form, 
and they form over a longer period of time, resulting in 
higher statistical fluctuations. 

A resolution study is shown in Figure El that compares 
the global Schmidt law computed with different numer- 
ical resolutions. Runs with different resolution converge 
within 10% in both the Ssfr and S gas . Although nu- 
merical resolution affects the total mass collapsed, and 
the number and location of fragments, as shown in Pa- 
per I, the SFR at tsf seems to be less sensitive to the 
numerical resolution. 

4.1. A Parameter Study 

In order to test how sensitive the global Schmidt law 
is to the radius R and the time t chosen to measure it, 
we carry out a parameter study changing both R and t 
individually. To maintain consistency with the previous 
section, we continue to assume a constant local SFE ee = 
30%. Figurc0compares the global Schmidt laws in total 
gas at different radii for the star-forming region R — 
i?50 and i?ioo (encircling 50% and 100% of the newly 
formed star clusters), while the time is fixed to t = tsf- 
We can see that the case with R 50 has larger scatter 
than that with i?ioo- This is due to the larger statistical 
fluctuations caused by the smaller number of star clusters 
within this radius. The global Schmidt law with R = 
.Rioo is almost identical to that with R = Rgo shown in 
Figure 

Figure |H1 compares the global Schmidt laws in total gas 
at different times t = 0.5tsf and t — 1.5tsf with the 
star formation radius fixed to R$q. Compared to the 
t = tsf case, models in the t = 0.5tsf case have higher 
Ssfr, because the SFR drops almost exponentially with 
time (§ 13. 1|) . Similarly, models in the t = 1.5tsf case 
shift to the lower right. Nevertheless, data derived from 
different times appear to preserve the power-law index of 
the Schmidt law, just differing in the n ormalizations. 

The global Schmidt law presented bv lLi et a (|2005b) 
was derived at a time when the total mass of the star 
clusters reached 70% of the maximum collapsed mass, 
which is close to 1.0 tsf in many models. The time in- 
terval At used to calculate the SFR was the time taken to 
grow from 30% to 70% of the maximum collapsed mass, 
rather than the At = 50 Myr used here. Nevertheless, 
the resu l ts pres ented here also agree well with those in 
ILi et al.l (|2005hlL 

Our parameter study demonstrates that the global 
Schmidt law depends only weakly on t he details of how it 
is measured. The small scatter seen in lKennicutti (0.998b) 
does suggest that additional physics not included in our 
modeling may be important. We should keep in mind 
that since we do not treat gas recycling, our models are 
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Fig. 5. — A comparison of the global Schmidt laws between our simulations and the obser vations. The red li ne is the least-square fit 
to the total gas of the simulated models, the black solid line is the best fit of observations from Kcnnicutt ( 1998b), while the black dotted 
lines indicate the observational uncertainty. The color of the symbol indicates the rotational velocity for each model (see Table 1); labels 
from M-l to M-4 are sub-models with increasing gas fraction; and open and filled symbols represent low and high T models, respectively. 



valid only within one gas consumption time tsf- The 
evolution after that may become unrealistic as most of 
the gas is locked up in the sinks. Nevertheless, our results 
suggest that the Schmidt law is a universal description 



of gravitational collapse in galactic disks. 

4.2. Alternative Global Star Formation Laws 
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Log E tot (M pc" 2 ) 

Fig. 6. — Same as Figure|2]but for the resolution study of low T 
models of G100-1 (blue) and G220-1 (orange). Models with total 
particle number of TVtot = 10 5 (Rl; open triangle), 8 X 10 5 (R8; 
open circle) and 6.4 X 10 6 (R64; open sguare) are shown. Models 
with regular resolution Ntot = 10 6 (RIO filled circle) are also shown 
for comparison. 



The existence of a well-denned global Schmidt law 
suggests that the star formation rate depends primar- 
ily on the g as surface density. As shown by several au- 
thors (e.g iQuirkl Il972t lLarsorJll98ct IKennicuttl Il989t 
lElmeereenl 119941 IKennicuttl Il998bj) . a simple picture 
of gravitational collapse on a free-fall timescale Tff oc 
p- 1 / 2 qualitatively produces the Schmidt law. Assum- 
ing the gas surface density is directly proportional to 
the midplane density, E gas oc p, it follows that Ssfr. oc 

Sgas/^ff Sgas- This suggests that the Schmidt law 
reflects the global growth rate of gas density under grav- 
itational perturbations. 

An alternative scenario that uses the local dynami- 
cal timescal e has been suggested by several groups (e.g. 
Wv^l 11983 IWvse fc Silldll987t ISilk] li"997t Klmegreli] 
19971 IHunter et aU 1199a iTanl 120001) In particular, 
Elmeereenl (Il997f ) and iHunter et all fl!998|) proposed a 
kinematic law that accounts for the stabilizing effect of 
rotational shear, in which the global SFR scales with the 
angular velocity of the disk, 

S , 

Ssfr oc — ^ oc E gas f2 (11) 

torb 

where t rb is the local orbital times cale and f2 is the or- 
bital frequency. IKennicuttl (|1998bf) gave a simple form 
for the kinematical law, 

Ssfr - 0.017 S gas O (12) 

with the normalization corresponding to a SFR of 21% 
of the gas mass per orbit at the outer edge of the disk. 




Log E tot (M pc 2 ) 

Fig. 7. — Same as Figure 151 but with different radii assumed 
for the star-forming region: (a) R = -R50, and (b) R = ftioo- The 
time is fixed at t = tsf- 



For our analysis, we follow IKennicuttl l)1998bft and de- 
fine t olh = 2ttR/V(R) = 2ir/Ct(R), where V(R) is the 
rotational velocity at radius R. We use the initial rota- 
tional velocity, which should not change much with time 
as it depends largely on the potential of the dark matter 
halo. Figure shows the relationship between Ssfr and 
Egas^ in our models. The densities of SFR Ssfr and to- 
tal gas S gas are calculated the same way as in Figure[S]at 
1.0 tsf and Rso, and fl(R) is calculated by using the ini- 



11 




-4 I , i , i , i , 

12 3 4 

Log S tot (M pc -2 ) 




Log E tot (M pc" 2 ) 

Fig. 8. — Same as Figure l5l but derived at different times (a) 
t = 0.5 Tgp and (b) t = f .5 rgp. The radius for the star formation 
region is fixed at i?80- 



tial total rotational velocity at Rsq . A least-square fit to 
the data gives E SFR = (0.036 ± 0.004) x (E gas n) 1 - 49±0 - 1 . 
This correlation has a steeper slope than the linear law 
given in equation l|12l) , suggesting a discrepancy between 
the behavior of our models and the observed kinematical 

law. 

1 I '■■> ver, iBoissier et al.l l|2003j) recently reported a 
slope of ~ 1.5 for the kinematical law from observations 



of 16 normal disk galaxies, in agreement with our results. 
Examination of Figure 7 in lKennicuttJ l|1998U) also shows 
that the normal galaxies, considered alone, seem to have 
a steepe r slope than the gal actic nuclei and starburst 
galaxies. Bois sier et al.l l)2003l) suggested several reasons 
for the discrepancy, the most important one being the 
difference bet ween their sample o f normal galaxies and 
the sample of [Kcnnicutt ( 1998b) including many star- 
burst galaxies and galaxy nuclei. More simulations, and 
models with higher SFR such as galaxy mergers are nec- 
essary to test this hypothesis. 

4.3. A New Parameterization 

The global Schmidt law describes global collapse in the 
gas disk. It does not seem to depend on the local star 
formation process. From Paper I and iLi et alJ (|2005b) 
we know that Ssfr correlates tight ly with the strength 
of gravitational instability (see alsolMac Low fc Klessenl 
l200l iKlessen et alJl200(l iHeitsch et alJl200lH . Here we 
quantify this correlation, using the Toomre Q parameter 
to measure the strength of instability. 

Figure ITUl shows the correlation between Esfr and the 
local gravitational instability parameters Q m in- The pa- 
rameters Qmin are minimum values of the Toomre Q pa- 
rameters for gas Qg,min(t), and the combination of stars 
and gas Q S g,mm{t) at a given time t, respectively . To ob- 
tain th e Q parameters, we follow the approach of lRaffkovl 
l)2001|) . as described in equations (l)-(3) of Paper I. We 
divide the entire galaxy disk at time t into 40 annuli, cal- 
culate the Q parameters in each annulus, then take the 
minimum. In the plots, the time when Ssfr is computed 
is t = Tgp . This correlation does not change significantly 
with time, but the scatter becomes larger at later times 
because the disk becomes more clumpy, which makes the 
calculation of Q g , m m(t) more difficult (see below). 

There is substantial scatter in the plots, at least partly 
caused by the clumpy distribution of the gas. Equations 
(l)-(3) for the Q parameters in Paper I are derived for 
uniformly distributed gas, such as in our initial condi- 
tions. As the galaxies evolve, the gas forms filaments 
or spiral arms probably leading to the fluctuations seen. 
The least-square fits to the data shown in Figure ITUl give 

Ssfr = (0.013 ± 0.003 M Q yr" 1 kpc~ 2 ) 

X [Q^mintTSF)]- 1 - 5 ^ - 23 (13) 

Ssfr = (0.019 ± 0.005 M Q yr" 1 kpc~ 2 ) 

X [Q^nintTSF)]- 1 - 12 "- 21 (14) 

If we take a first-order approximation, Q sg oc 1 /E gas , 
then equation ifHfl) gives Esfr oc Eg a g 4 at t = rsp, agree- 
ing very well with the observations. The slopes derived 
from Q g appear to be lower than those derived from Q sg , 
but are still within the slope range observed. 

Keep in mind that the local instability is a non-linear 
interaction between the stars and gas, and so is much 
more complicated than the linear stability analysis pre- 
sented here. Also, the instability of the entire disk at a 
certain time is not fully represented by the minimum val- 
ues of the Q parameters we employ here, although they 
do represent the region of fastest star formation. These 
factors limit our ability to derive the global Schmidt law 
directly from the instability analysis. 
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5. LOCAL SCHMIDT LAW 

The relationship between surface density of SFR Ssfr 
and gas density S gas can also be measured as a function 
of radius within a galaxy, givi n g a lo cal S chmidt law. 
Observations by iWong Ez Blitd l|2002|) and iHever et al.1 



( 2004) show significant variations in both the indices N 
and normalizations A of the l ocal Schmidt l aws o f indi- 
vidual galaxies. For example, IHever et al.l l)2004|) show 
that M33 follows the law 
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J SFR 



M 



yr 



1 kpc- a ) 

. 3.3±0.07 



(0.0035 ±0.066 

x (E tot /l M pc-y™' , (15) 

while iWong fc Blitd H002) show that N has a range of 
1.23-2.06 in their sample. 

To derive local Schmidt laws, we again divide each indi- 
vidual galaxy into 40 radial annuli within 4 i?<j, then com- 
pute E gas and EgpR. in each annulus. The SFR is mea- 
sured ai t — rgp as in §0] For models where rgp > 3 Gyr 
or is beyond the simulation duration, the maximum sim- 



ulated timestep is used instead, as listed in Table 1. 
5.1. Star Formation Thresholds 

Figure ^] shows the relation between EgpR and E gas 
correlations for all models in the simulations that form 
stars in the first 3 Gyr. With such a large number of 
models in one plot, it is straightforward to character- 
ize the general features, and t o compare with th e obser- 
vations shown in Figure 3 of IKenni cuttJ (ll 998bD. Simi- 
lar to the individual galaxi es in IKennicuttJ l)1998b|) and 
iMartin fc KennicuttJ 1)200 ll) . each model here shows a 
tight EsFR-Eg as correlation, the local Schmidt law. How- 
ever, Esfr drops dramatically at some gas surface den- 
sity. This is a clear indication of a star-formation thresh- 
old. 

We therefore define a threshold radius i?th as the ra- 
dius that encircles 95% of the newly formed stars. The 
gas surface density at the threshold radius i?th in Fig- 
ure ^2 has a range from ~ 4M Q pc~ 2 for the relatively 
stable model G220-1 (low-T) to - 60M Q pc~ 2 for the 
most unstable model G220-4 (high-T) . Note that in some 
galaxies, there are also small e r dip s in SFR at higher 
density. IMartin fc KennicuttJ 1)20011) suggest that rota- 
tional shearing can cause an inner star-formation thresh- 
old. However, the inner dips in our simulations are likely 
due to the lack of accretion onto sink particles in the sim- 
ulations after most of the gas in the central region has 
been consumed. Further central star formation in real 
galaxies would occur due to gas recycling, which we ne- 
glect, and, probably more important, after interactions 
with other galaxies. 

In the analysis of observations, a dimensionless param- 
eter, aQ = E t h/E cr it = 1/Q has been introduced to re- 
late the star formation t hreshold to the Toomre unstable 
radius l)KennicuttJ 1 1 989|) . The critical radius is usually 
defined as the ra dius where Q„ = 1. With a sample of 
15 spiral galaxies. iKemiiciifrn 1)198 91 found aQ ~ 0.63 by 
assuming a constant effective sound speed (the velocity 
dispersion) of the gas c« = 6 km s -1 . T his result was 
confirmed by IMartin fc Kennicuttl l)2001l) with a larger 
sample of 32 well-studied nearby spiral galaxies, who re- 
ported a range of ap ~ 0.3- 1 .2, wi th a median value of 
0.69. However. iHunter et alJ l)1998j) found a Q ~ 0.25 for 
a sample of irr egular galaxie s with c s = 9 km s _1 . As 
pointed out bv lSchaveIl)2004l) . this derivation of aQ de- 
pends on the assumption of c s . The values of aQ derived 
from our models using their actual values of c s as shown 
in Figure El We find that the value of «q depends not 
only on the gas sound speed, but also on the gas fraction 
of the galaxy. For models with the same rotational ve- 
locity and gas fraction, lower gas sound speed results in a 
higher value of aQ. For models with the same total mass 
and sound speed, higher gas fraction leads to higher aQ . 
The gas-poor models in our simulations (/ g = 20%) have 
a range of aQ ~ 0.2-1.0, agreeing roughly with observa- 
tions. This again may reflect the relative stability of the 
nearby galaxies in the observed samples. 

There are several theoretical approaches to ex- 
plain the presence of s tar formation thresholds. 
IMartin fc Kennicuttl 1)200 U) suggest that the gravita- 
tional instability model explains the thresholds well, 
with the deviation of aQ from one simply due to 
the n on-uniform di s tribut ion of gas in real disk galax- 
ies. IHunter et alJ l)1998f) proposed a shear criterion 
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for star forming dwarf irregular galaxies, as t hey ap- 
pear t o be sub-critical to the Toomre criterion. iSchavd 
(2004) modeled the thermal and ionization structure of 
a gaseous disk. He found the critical density is about 
^crit ~ 3-10 Mq pc -2 with a gas velocity dispersion 



of ~ 10 km s _1 , and argued that thermal instability 
determines the star formation threshold in the outer 
disk. Our models suggest that the threshold depends 
on the gravitational instability of the disk. The de- 
rived S cr ; t and otQ from our stable models (Q S g,min > 1) 
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a gree well with observations, supporting the arguments 
of IMartin fe Kennicuttl l |200lj) . 

5.2. Local Correlations Between Gas and Star 
Formation Rate 

We fit the local Schmidt law to the total gas surface 
density within Rth, as demonstrated in Figure HHl The 
models in Figure EI] all have the same rotational veloc- 
ity of 220 km s _1 but different gas fractions and sound 
speeds. The local Schmidt laws of these models vary only 
slightly in slope and normalization. 

Figures 114b and 114b compare the slope N and nor- 
malization A of the local Schmidt laws for all models in 
Table 1 that form stars in the first 3 Gyr. The slope of 
the fit to the total gas in Figurefrik varies from about 1.2 
to 1.7. Larger galaxies tend to have larger N. However, 
the average slope is around 1.3, agreeing reasonably well 
with that of the global Schmidt law. 

There is substantial fluctuation in the normalization of 
fits to the total g shown in Figure IT4T &'). The vari- 
ation is more than an order of magnitude, with gas-rich 
models tending to have high A. However, the average 
value of A settles around 2.2, agreeing surprisingly well 
with that of the global Schmidt law. Overall, the aver- 
aged local Schmidt law gives: 

-4 



Ssfr = (2.46 ± 1.62 x 1CT 4 M yr" 1 kpc~ 2 

/ V \ 1-31±0.15 



v l M Q pc- 



(16) 



The averaged local Schmidt law is very close to the 
global Schmidt law in § 0] 



The average slope in equation 11611 is rather smaller 
than the value of N = 3.3 observed bv lHever et alJ l|2004D 
in M33. The galaxy M33 is very interesting. It is a nearly 
isolated, small disk galaxy with low luminosity and low 
mass. It has total mass of M to t ~ 10 11 and a gas 
mass of M gas - 8.0 x 10 9 M Q l|Hever et all 12004]) . It 
is molecule-poor and sub-critical, with gas surface den- 
sity is much smaller than t he threshold surface densit y 
for star formation found bv IMartin fe Kennicuttl 1 2001). 
However, it is actively forming stars ()Heve^et*^Ll] 200^) . 
We do not have a model that exactly resembles M33, al- 
though a close one might be model G 100-1 in terms of 
mass. However, the gas velocity dispersion of M33 is un- 
known, so we cannot make a direct comparison with our 
G100-1 models. In Figured the low-T model G100-1 
has N ~ 1.4, but we have not derived a value for its 
high-T counterpart, as it does not form stars at all in 
the first 3 Gyrs. Any stars that form in a disk similar 
to this will likely form in spiral arms or other nonlinear 
density perturbations that are not well characterized by 
an azimuthally averaged stability analysis. If these per- 
turbations occur in the highest surface density regions 
as might be expected, the local Schmidt law will have a 
very high slope as observed. This speculation will need 
to be confirmed with models reaching higher mass reso- 
lution in the future. The details of the feedback model 
and equation of state may also begin to play a role in 
this extreme case. 

The averaged values of our derived loca l Schmidt laws 
do ag ree well with the observations by iWong fc Blitzl 
( 2002) of a number of other nearby galaxies. The similar- 
ity between the global and local Schmidt laws suggests a 
common origin of the correlation between Ssfr and S gas 
in gravitational instability. 

6. STAR FORMATION EFFICIENCY 

The SFE is poorly understood, because it is difficult 
in both observations and simulations to determine the 
timescale for gas removal and the gaseous and stellar 
mass within the star formation region. On the molec- 
ular cloud scale, observations of several nearby embed- 
ded clusters with mass M < 1000M Q ind icate that the 
SFEs range from approximately 10-30% l)Lada fc Ladal 
2003). However, it is thought that field stars form 
with SFE of only 1-5% in gi ant molecular clouds (e.g., 
iDuerr. Imhoff. fc Ladal Il982|) . while the formation of 
a bound stellar cluster requires a local SFE > 20- 
50% (e.g. TWilking fc Lac^ll983t lEhnegreen fc Efremovl 
19971). An analytical model including outflows by 
Matzner fc McKeeT (|2000^ suggests that the efficiency of 
cluster formation is in the range of 30-50%, and that 
of single star formation could be anywhere in the range 
25-70%. 

In the analysis of our simulations presented here, we 
convert the mass of the sink particles into stars using a 
fixed local SFE ei = 30%, consistent with both the ob- 
servations and theoretical predictions mentioned above. 
This local efficiency is different from the global star for- 
mation efficiency in galaxies e g < ei, which measures 
the fraction of the total gas turned into stars. On a 
galactic scale, the star formation efficiency appears to 
be associated with th e fraction of molecular gas (e.g., 
iRownd fc YoTmlll999T) . The global SFE has had values 
derived from observations over a wide range, depend- 
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ing o n the gas distribu tion and the molecular gas frac- 
tion l)KennicuttJll998bT) . For example, in normal galaxies 
e g ~ 2-10%, while in starburst galaxies e g = 10-50%, 
with a median value of 30%. One factor that appears 
to contribute to the differences in e g is the gas content. 



The global SFE is generally averaged over all gas compo- 
nents, but since star formation correlates tightly with the 
local gravitational instability one expects higher global 
SFE in more unstabl e galaxies. In fact, as pointed out by 
IWong fc Blitd l|2002[K most normal galaxies in the sample 



17 



1.8 



1.6 



1.4 



1.2 



1.0 



100.0 



10.0 



1.0 



0.1 



□ 


LT:M 




A 


LT:M 


-2 


<> 


LT:M 


-3 


O 


LT:M 


-4 



■ 


HT 


M-1 


A 


HT 


M-2 


♦ 


HT 


M-3 


• 


HT 


M-4 



o 

A 



♦ 



50 100 150 200 250 
V 200 (km s" 1 ) 



-i — i — i — r- 



- 1 — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i - 



□ 


LT:M 




A 


LT:M 


-2 





LT:M 


-3 


o 


LT:M 


-4 



■ 


HT:M 




A 


HT:M 


-2 


♦ 


HT:M 


-3 


• 


HT:M 


-A 



o 



o 

A 



♦ • 

♦ 

□ 



♦ 

t 



_i_ 



50 100 



'200 



(km 



150 



200 250 



Fig. 14. — Local Schmidt law (a) slope N and (6) normalization 
A for all models with r$p < 3 Gyr. The fit is to the total gas 
surface density Stot within the threshold radius R t ^ . The symbols 
represent different models, as indicated in the legends. The dotted 
line in each panel indicates the linearly averaged value across the 
models shown. 



of lKennicutd Jl998b) are molecule-poor galaxies, which 
seem to have high stability and low SFE, while molecule- 
rich starburst galaxies appear to be unstable, forming 
stars with high efficiency. 

The variation of the normalization of the local Schmidt 
laws in § also suggests that the global SFE varies from 
galaxy to galaxy. To quantitatively measure the SFE in 
our models, we apply the common definition of the global 
SFE, 

e g = M* / (M* + A/gas ) = M* / M i„it (17) 



over a period of 10 8 years, an average timescale for star 
formation in galaxy. In this equation, M* = QM s ; n k is 
the mass of newly formed stars, and M gas includes both 
the remaining mass of the sink particles and the SPH 
particles, so that M* + M gas = M; n i t is the total mass of 
the initial gas. 

Figure I15f a) and ( b) show the relation between the 
minimum values of the initial Q parameter Q 



s<?,min 



and 

the global SFE normalized by the local SFE e g /ee. The 
time period is taken as the first 100 Myr after star for- 
mation starts. If we take et as a constant for all models, 
it appears that e g declines as Q S p,min increases. There- 
fore, e g is high in less stable galaxies with high mass or 
high gas fraction. A least-absolute-deviation fit of the 
data gives a linear fit of e g /e£ = 0.9 — 0.97 Q S g,min- This 

< 1. For more stable 
the SFE remains 



fit is good for values of Q 



S3, mm 



galaxies, with larger values of Q S g,mm, 
finite, deviating from the linear fit. 

Using the empirical relations we have derived from our 
models earlier in the paper, we can derive a better ana- 
lytic expression for e g . Equation i|17|) can be combined 
with equation J3Jl in § 13.11 to yield an equation for the 
global SFE 



e g = M [1 - cxp (-*/t S f)] 



(18) 



We evaluate this at t = 100 Myr, taking the defini- 
tions of Mo and tsf derived from equations (@J and (JjJJ. 
Normalizing by the local SFE, we find 



,/e e = 0.96 x [1 - 2.88exp(-1.7/Q sg , mi „)] 
x [l - exp (-2.9e- Q "/ ' 24 



(19) 



The function given by equation (|19fl is shown in Fig- 
ure El For Qs^min < 1.0 it is well approximated by the 
much simpler linear function 

e g /e e ~ 0.9 - Qsg^n (20) 

as shown in Figure irKT a - ). At larger values of Qs^mm, the 
exact function predicts the SFE in our models excellently 
as shown in Figure H5T M. Observational verification of 
this behavior is vital. 

7. ASSUMPTIONS AND LIMITATIONS OF THE MODELS 

7.1. Isothermal Equation of State 

One of the two central assumptions in our model is 
the use of an isothermal equation of state to represent 
a constant velocity dispersion. This is of course a sim- 
plification, as the interstellar medium in reality has a 
broad range of temperatures 10 K< T < 10 7 K. However, 



neutral gas velocity dispersions in normal spi ral galaxies 
:over a f ar more limited range, as reviewed bvlKennic 
<1998rj|) : lETmegreen fc Scalo I pOO^MScalo k Elmegr 



cover a f ar more limited range, as reviewed byl Kennicutt 
, _ , _ jcc:i 

(|2004f) and lDib^t^dTl|2005li . The characteristic a in- 
creases when the averaged Ssfr of a galaxy reaches 
tens of solar masses per year, but the normal galax- 
ies in the sample with re liable measurements lie in the 
range ~7-13 k m s" 1 fe.g..lElmegreen fc Elmegreer~ 
i Meurer et all 119961: Ivan Zee et al.l 119971 iStil fc Israeli 
12002 IHinnelein et al. 1120031) . 

At least two mechanisms appear viable for main- 
taining roughly constant velocity dispersion for the 
bulk of the gas in a galactic disk, supernova feedback 
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and magnetorotational instability l)Mac Low fe Klessenl 
I2004|) . Three-dimensional simulations in a periodic box 
with para meters characteris tic of the outer parts of galac- 
tic disks by Di b et all ((2005) show that supernova driving 
leads to constant velocity dispersions of a ~ 6 km s _1 
for the total gas and cthi ~ 3 km s _1 for the Hi 
gas, independent of the supernova rate. Simulations of 
the feedback effects across whole galactic disks do sug- 
gest that the inn er parts have slightly higher velocity 
dispersions (e.g. iThacker fc Couchmanl 12000). though 
within the range that we consider. The magnetoro- 



tational instability in ga lactic disks was suggested by 
iSellwood fc Balbiia l)1999ft to maintain the observed ve- 
locity dispersion, a sugge stion that has since been sub- 
stantia ted by both local fPiontek fc Ostrike'rll2005tl and 
global (Dziourkcvitch. Elstncr. & Riidiger 2004) numer- 
ical models. This may act even in regions with little or 
no active st ar formation. 

Recently, iRobertson et al.l l)2004|) presented simula- 
tions of galactic disks and claimed that an isothermal 
equation of state leads to a collapsed disk as the gas 
fragments into clumps that fall to the galactic center 
due to dynamical friction. Howe ver, similar behavior is 
seen in models bv llmmeli et al. 1 1 



( 200J) who did not use 
an isothermal equation of state, but al so ran at resolu- 
tions not satisfying the Jea ns criterion l)Bate fc BurkertJ 
119971 iTruelove etafl I1997|) . On the other hand, us- 
ing essentially t he sa me code and galaxy model as 
IRobertson et alJ ((2004) , but with higher resolution sat- 
isfying the Jeans criterion, we do not see this collapse. 
Insufficient resolution that fails to resolve the Jeans mass 
leads to spurious, artificial fragmentation and thus col- 
lapse. 

Sin : atkn by IGovernato et al.l l)2004j) suggest that 
some long-standing problems in galaxy formation such 
as the compact disk and lack of angular momentum may 
well be due to insufficient resolution or violation of nu- 
merical criteria. Our results lead us to agree that the 
isothermal equation of state is not the cause of the com- 
pact disk problem, but rather inadequate numerical res- 
olution. 

Our assumption of an isothermal equation of state 
does, of course, rule out the treatment by our model 
of phenomena such as galactic winds associated with the 
hot phase of the interstellar medium (although the vent- 
ing of supernova energy vertically may help maintain the 
isothermal behavior of the gas in the plane) . The strong 
starbursts produced in some of our galaxy models will 
certainly cause strong galactic winds. It remains unclear 
whether even strong starbursts can remove substantial 
amounts of gas, though. Certainly they cannot in small 
galaxies (|Mac Low fc Ferraralll999j) . and larger galaxies 
would seem more resistant to stripping in starbursts than 
smaller ones. However, galactic winds will certainly in- 
fluence the surroundings of starburst galaxies, as well as 
their observable properties. These effects should eventu- 
ally be addressed in future simulations with more com- 
prehensive gas physics and a more realistic description of 
the feedback from star formation. 

7.2. Sink Particles 

The use of sink particles enables us to directly identify 
high gas density regions, measure gravitational collapse, 
and follow the dynamical evolution of the system to a 
long time. We can therefore determine the star formation 
morphologies and rates, and study the Schmidt laws and 
star formation thresholds. 

However, one shortcoming of our sink particle imple- 
mentation is that we do not include gas recycling. Once 
the gas collapses into the sinks, it remains locked up 
there. As discussed in Paper I, the bulk of the gas that 
does not form stars will remain in the disk and contribute 
to the next cycle of star formation. Also, the ejected ma- 
terial from massive stars will return into the gas reservoir 
for future star formation. Another problem is accretion. 
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In the current model, sink particles accrete until the sur- 
rounding gas is completely consumed. However, in real 
star clusters, the accretion would be cut-off due to stellar 
radiation, and the clusters will actually lose mass due to 
outflow and tidal stripping. 

These shortcomings of our sink particle technique may 
contribute to two limitations of our models: first, the 
decline of star formation rate over time due to consump- 
tion of the gas, as seen in Figure^ second, the variation 
of SFR over time in the simulated global Schmidt laws 
in Figure |S| Nevertheless, as we have demonstrated in 
the previous sections, our models are valid within one 
gas consumption time rsp, and are sufficient to inves- 
tigate the dominant physics that controls gravitational 
collapse and star formation within that period. 

7.3. Initial Conditions of Galaxies 

Many nearby galaxies appear to be gas-poor and stable 
(Qsg > 1)- However, their progenitors at high redshift 
were gas rich, so the bulk of star formation should have 
taken place early on. In order to test this, we vary the gas 
fraction (in terms of total disk mass) in the models. We 
also vary the total galaxy mass and thus the rotational 
velocity. These result in different initial stability curves. 
Massive, or gas-rich galaxies have low values of the Q 
parameters, so they are unstable, forming stars quickly 
and efficiently. 

There are no observations yet that directly measure 
the Q values in starburst galaxies. How ever, indirectly, 
observations bv iDalcanton et aD l)2004l) show that dust 
lanes, which trace star formation, only form in unsta- 
ble regions. Mo reover, observations of c olor gradients in 
disk galaxies bv lMacArthur et alJ l)2004f) show that mas- 
sive galaxies form stars earlier and with higher efficiency. 
Both of these observations are naturally explained by our 
models. 

The Toomre Q parameter for gas Q g differs from that 
for a combination of stars and gas Q sg in some of our 
model galaxies. This leads to slightly different results in 
Figure |3] and Figure ^| where we compare Q S g,mm and 
Qg,min- However, we believe Q sg is a better measure 
of gravitational instability in the disk, as it takes into 
account both the collisionless and the collisional compo- 
nents, and the interaction between them. We note that it 
is a simplified approach to quantify the instability of the 
entire disk with just a number Q S g )m i n , as Q has a radial 
distribution, evolves with time, and is an azimuthally 
averaged quantity, but nevertheless, we find interesting 
regularities by making this approximation. 

8. SUMMARY 

We have simulated gravitational instability in galaxies 
with sufficient resolution to resolve collapse to molecu- 
lar cloud pressures in models of a wide range of disk 
galaxies with different total mass, gas fraction, and ini- 
tial gravitational instability. Our calculations are based 
on two approximations: the gas of the galactic disk has 
an isothermal equation of state, representing a roughly 
constant gas velocity dispersion; and sink particles are 
used to follow gravitationally collapsed gas, which we as- 
sume to form both stars and molecular gas. With these 
approximations, we have derived star formation histo- 
ries; radial profiles of the surface density of molecular and 
atomic gas and SFR; both the global and local Schmidt 



laws for star formation in galaxies; and the star forma- 
tion efficiency. 

The star formation histories of our models show the 
exponential dependence on time given by equation |6] in 
agreement with, for e xample, the interpretatio n of galac- 
tic color gradients bv iMacArthur et alJ l)2004l) . The ra- 
dial profiles of atomic and molecular gas qualitatively 
agree with those observed in nearby galaxies, with sur- 
face density of molecular gas peaking cen trally at values 
much above that of the atomic gas (e.g. iWong fc Blitd 
2002). The radial profile of the surface density of SFR 
correlates linearly with that of the molecular gas, ag ree- 
ing with the observations of iGao fc Solomon! (|2004bT l. 

Our models quan titatively reprodu ce the observed 
global Schmidt law <)Kennicuttill998bj) — the correlation 
between the surface density of star formation rate Ssfr 
and the gas surface density S gas — in both the slope and 
normalization over a wide range of gas surface densi- 
ties (eq. EJ)- We show that Esfr is strongly correlated 
with the gravitational instability of galaxies Esfr cx 

[Q S9 ,min(TsF)] _1 ' 54±0 ' 23 , where Q S9 , m in(TSF) is the local 
instability parameter at time t = rsp (see eq. II 31) . This 
correlation naturally leads to the Schmidt law. 

On the other hand, our models do not reproduce the 
corr elation Ssfr ~ Sg as Q derived from kinematical mod- 
els l)Kennicuttl 1 1 998bF . However, they may agree better 
with the dependence of the normal galaxie s on this quan- 
tity, as suggested bv iBoissier et all (|2003Tl . The discrep- 
ancy may be caused by the lack of extreme starburst 
galaxies such as galaxy mergers in our set of models. 

The local Schmidt laws of individual galaxies clearly 
show evidence of star formation thresholds above a crit- 
ical surface density. The threshold surface density varies 
with galaxy, and appears to be determined by the grav- 
itational stability of the disk. The derived threshold 
parameters for our stable models cover the range of 
values seen in observations of normal galaxies. The 
local Schmidt laws have significant variations in both 
slope and normaliz ation, but also cover the observational 
ranges rep orted bvlWong fc Blitz! l|2002fi . IBoissier et alJ 
pOffo and IHever et a,U <|2004D . The average normaliza- 
tion and slope of the local power-laws are very close to 
those of the global Schmidt law. 

Our models show that the global star formation effi- 
ciency (SFE) t g can be quantitatively predicted by the 
gravitational instability of the disk. We have used a fixed 
local SFE q = 30% to convert the mass of the sink par- 
ticles to stars in our analysis. This is a reasonable as- 
sumption for the SFE in dense, high pressure molecular 
clouds. The global SFE of a galaxy then can be shown 
to depend quantitatively on a nonlinear function (eq. I19f) 
of the minimum Toomre parameter Q S g,min for stars and 
gas that can be approximated for Q S g,min < 1.0 with the 
linear correlation e g /ei cx 0.9 — Q S g,min- More unstable 
galaxies have higher SFE. Massive, or gas-rich galaxies 
in our suite of models are unstable, forming stars quickly 
with high efficiency. They represent starburst galaxies. 
Small, or gas-poor galaxies are rather stable, forming 
stars slowly with low efficiency, corresponding to quies- 
cent, normal galaxies. 
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TABLE 1 

Galaxy Models and Numerical Parameters 



Model a 


/g b 




<9s 9 (LT) d 


Q S9 (HT) C 


S T f 

AW 




m g 


tsfCLT) 1 


T SF (HT)j 


G50-1 


0.2 


1.41 


1.22 


1.45 


1.0 


10 


0.08 


4.59 




G50-2 


0.5 


1.41 


0.94 


1.53 


1.0 


10 


0.21 


1.28 




G50-3 


0.9 


1.41 


0.65 


1.52 


1.0 


10 


0.37 


0.45 




G50-4 


0.9 


1.07 


0.33 


0.82 


1.0 


10 


0.75 


0.15 


0.53 


G100-1 


0.2 


2.81 


1.08 


1.27 


6.4 


7 


0.10 


2.66 




G100-2 


0.5 


2.81 




1.07 


1.0 


10 


1.65 






G100-3 


0.9 


2.81 




0.82 


1.0 


10 


2.97 




1.92 


G100-4 


0.9 


2.14 




0.42 


1.0 


20 


5.94 




0.15 


G120-3 


0.9 


3.38 




0.68 


1.0 


20 


5.17 




0.46 


G120-4 


0.9 


2.57 




0.35 


1.0 


30 


10.3 




0.16 


G 160-1 


0.2 


4.51 




1.34 


1.0 


20 


2.72 




3.1 k 


G160-2 


0.5 


4.51 




0.89 


1.0 


20 


6.80 




0.58 


G160-3 


0.9 


4.51 




0.52 


1.0 


30 


12.2 




0.30 


G160-4 


0.9 


3.42 




0.26 


1.5 


40 


16.3 




0.11 


G220-1 


0.2 


6.20 


0.65 


1.11 


6.4 


15 


1.11 


0.28 


3.0 k 


G220-2 


0.5 


6.20 




0.66 


1.2 


30 


14.8 




0.39 


G220-3 


0.9 


6.20 




0.38 


2.0 


40 


15.9 




0.25 


G220-4 


0.9 


4.71 




0.19 


4.0 


40 


16.0 




0.096 



a First number is rotational velocity in km s _1 at the virial radius, the second number indicates 
sub-model. Sub-models have varying fractions mj of total halo mass in their disks, and given 
values of / g . Sub-models 1-3 have = 0.05, while sub-model 4 has m& = 0.1. 

b Fraction of disk mass in gas. 

c Stellar disk radial exponential scale length in kpc 

d Minimum initial value of Q sg {R) for low-T models. 

c Minimum initial value of Q ag (R) for high-T models. 

'Total particle number in units of 10 6 

g Gravitational softening length of gas in pc. 

h Gas particle mass in units of 10 4 Mq. 

'Star formation timcscale in Gyr of low-T model (from Paper I). 
J Star formation timescale in Gyr of high-T model (from Paper I). 
k Maximum simulation timestep instead of the star formation timescale tsf- 



